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Abstract 

The results obtained by several experiments on atmospheric neutrino oscillations are summarized and 
discussed. Then the results obtained by different long baseline neutrino experiments are considered. Finally 
conclusions and perspectives are made. 

1 Introduction 

Atmospheric neutrinos are well suited for the study of neutrino oscillations, since they have energies from a 
fraction of GeV up to more than 100 GeV and they may travel distances L from few tens of km up to 13000 
km; thus LjE v ranges from ~1 km/GeV to 10 5 km/GeV. Atmospheric neutrinos may study v oscillations 
for small Am 2 . 

The early water Cherenkov detectors 1MB [1] and Kamiokande [2] reported anomalies in the ratio of 
muon to electron neutrinos, while the tracking calorimeters NUSEX [3] and Frejus [4], and the Baksan 
[5] scintillator detector did not find any. In 1995 MACRO found a deficit for upthrougoing muons [6]. 
Then the Soudan 2 experiment [7] confirmed the ratio anomaly. In 1998 Soudan 2 [8], MACRO [9] and 
SuperKamiokande (SK) [10] reported deficits in the fluxes with respect to Monte Carlo (MC) predictions 
and angular distribution distortions; instead the v e distributions were in agreement with non oscillated 
MCs. These features may be explained in terms of v ^ < — ► v T oscillations. 

The atmospheric neutrino flux was computed by many authors in the mid 1990s [11] and in the early 
2000s [12]. The last had improvements, but also a new scale uncertainty. 

Several long baseline v beams were and are operational: KEK to Kamioka (K2K) (250 km), NuMi from 
Fermilab to the Soudan mine (730 km) and CERN to Gran Sasso(GS) (CNGS) (730 km). The K2K [13] 
and MINOS [14] experiments obtained results in agreement with the atmospheric v results. The CNGS 
beam is under tests: some events were recorded in OPERA [15] and LVD [16] at GS. 

Most experiments are disappearance experiments; OPERA is an appearance experiment, searching for 
v T in a pure v ^ beam. 

Long baseline experiments will give an increasing contribution to neutrino physics. 



2 Atmospheric neutrino oscillations 

A high energy primary cosmic ray (CR), proton or nucleus, interacts in the upper atmosphere producing 
a large number of charged pions and kaons, which decay yielding muons and z^'s; also the muons decay 
yielding and v e . The ratio of the numbers of to v e is ~ 2 and N v jNy ~ 1. Atmospheric neutrinos are 
produced at 10-20 km above ground, and they proceed towards the earth. 

If z/s have non-zero masses, one considers the weak flavour eigenstates v e , v ^ , v T and the mass 
eigenstates v\, ^2, ^3. Flavour eigenstates are linear combinations of mass eigenstates. For 2 flavour 
(^/i, v T ) and 2 mass eigenstates (zv 2 , ^3) one writes 

f = V2 COS #23 + ^3 Sin #23 

\v T = -v 2 sin #23 + ^3 cos #23 
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Figure 1: a) Cross section of the MACRO detector, sketch of event topologies, b) Schematic layout of the 
SuperKamiokande detector 



where #23 is the mixing angle. The survival probability of the "beam" is 

/ 1 07 Am 2 . T \ 

P("» - = 1 - - v T ) = 1 - sin 2 20 2 3 sin 2 ( — J (2) 

where Am 2 = m\ — m\ and L is the distance from v production to detection. 

The simple formula Eq. 2 is modified by additional flavours and by matter effects. 

3 Results from atmospheric neutrino experiments 

After the early indications of an anomaly in atmospheric z/s [1-7], in 1998 Soudan 2, MACRO and SK 
provided strong indications in favour of < — > v T oscillations [8-10] [17]. Confirming results were presented 
by long baseline experiments [13] [14]. 

Soudan 2 used a modular fine grained tracking and showering calorimeter of 963 t located underground 
in the Soudan Gold mine in Minnesota. The detector was made of lmxlmx2.5m modules, surrounded by 
an anticoincidence shield. The bulk of the mass was 1.6 mm thick corrugated steel sheets interleaved 
with drift tubes [8]. The Soudan 2 double ratio for the zenith angle range —1 < cos 6 < 1 is 
R' = (N^ I N e ) data I (Nft I N e ) mc = 0.68 ± 0.ll stat , consistent with oscillations. 

MACRO detected upgoing v^s via CC interactions — > \i\ upgoing muons were identified with 
streamer tubes (for tracking) and scintillators (for time-of-flight measurements), points and rectangles in 
Fig. la [18]. Events were classified as: 

Upthroughgoing muons from interactions in the rock below the detector of with (E u ) ~ 50 GeV. Data 
were compared with the predictions of the Bartol96 [11], FLUKA, HKKM01 [12] MCs, Fig. 2a. The 
shape of the angular distribution and the absolute value favoured v oscillations with Am 2 = 0.0025 eV 2 
and maximum mixing. The absolute value of the muons is 25% higher than FLUKA and HKKM01 MC 
predictions . The difference between new and old MC predictions is due to new fits of CR data [21]. 

Low energy events. Semicontained upgoing muons (IU) come from interactions inside the lower 
apparatus. Up stopping muons (UGS) are due to external interactions yielding upgoing muons stopping 
in the detector; the semicontained downgoing muons (ID) are due to downgoing v^s interacting in the lower 
detector; the lack of time information prevents to distinguish between the two samples. An equal number 
of UGS and ID events is expected. The average parent neutrino energy is 2-3 GeV. The data are compared 
with MC predictions without oscillations in Figs. 2b, c: they show a uniform deficit over the whole angular 
distribution with respect to Bartol96 predictions. 
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Figure 2: left) MACRO upthroughgoing muons compared with the oscillated MC predictions of Bartol96 
(solid curve), HKKM01 (dash-dotted line), FLUKA fitted to the new CR fit (dashed curve) and FLUKA 
with the old CR fit (dotted curve), center) IU and right) ID+UGS events (black points) compared with the 
non oscillated Bartol96 MC (dashed line) and < — ► v T predictions. 



< — > v T against < — > v 'sterile- Matter effects due to the difference between the weak interaction 
effective potential for with respect to sterile neutrinos, yield different total number and different zenith 
distributions of upthroughgoing muons. R m eas between the events with —1 < cosO < —0.7 and with 
—0.4 < cos 9 < was used [9]. The measured ratio, Rmeas = 1-38, is compared with R T = 1.61 and 
Rsterile = 2.03. One concludes that < — > v steri i e oscillations are excluded at the 99.8% c.l. 

v^x energy estimate by Multiple Coulomb Scattering (MCS) of upthroughgoing muons. An estimate of 
the muon energy was made through their Multiple Coulomb Scattering (MCS) in the absorbers [19]. The 
ratios Data/MCno osc as a function of (L/E u ) are in agreement with < — > v T oscillations [20]. 

New determination of the oscillation parameters. In order to reduce the effects of systematic uncertainties 
in the MCs, MACRO used the following three independent ratios (it was checked that FLUKA, HKKM01 
and Bartol96 MC yield the same predictions): 

(i) High Energy Data: zenith distribution ratio: R\ = N vert /Nh or 

(ii) High Energy Data: energy measurement ratio: R2 = Ni ow /Nhi g h 
(hi) Low Energy Data: R 3 = (Data/ MC) W / (Data/ MC) ID+UGS . 

The no oscillation hypothesis had a probability P~3 -10 -7 and is ruled out by ~ 5a. Fitting the 3 ratios to 
the v ^ < — ► v T oscillation formulae, one obtained sin 2 26 = 1, Am 2 = 2.3 • 10 -3 eV 2 . Using Bartol96, one 
adds the information on absolute fluxes: 

(iv) High energy data (systematic error ~17%): Ra = N meas /NMC' 

(v) Low energy semicontained muons (scale error 21%): R§ = N meas /NMC' 
These leave the best fit values unchanged and improve the significance to 6a. 

SuperKamiokande (SK) is a large cylindrical water Cherenkov detector of 39 m diameter and 41 m 
height containing 50 kt of water (fiducial mass 22.5 kt); it is seen by 50-cm-diameter inner-facing phototubes 
(PMTs), Fig. lb. The 2 m thick outer layer of water acts as an anticoincidence; it is seen by smaller outward- 
facing PMTs. The detector is located in the Kamioka mine, Japan. Atmospheric neutrinos are detected 
in SK by measuring the Cherenkov light generated by the charged particles produced in the neutrino CC 
interactions with water nuclei. Thanks to the high PMT coverage, the experiment detects events of energies 
as low as ~5 MeV. 

The large detector mass allows to collect a high statistics sample of fully contained events (FC) up to 
~5 GeV. The FC events yield rings of Cherenkov light on the PMTs. FC events can be subdivided into 
sub-GeV and multi-GeV events, with energies below and above 1.33 GeV. events include only single-ring 
events, while multi-ring ones (MRING) are treated as a separate category. The partially contained events 
(PC) are CC interactions with vertex within the fiducial volume and at least a charged particle, typically 
the /i, exits the detector (the light pattern is a filled circle). Upward-going muons ( UPMU), produced by 
from below interacting in the rock, are subdivided into stopping ((E u ) ~ 7 GeV) and throughgoing muons 
((E v ) - 70^80 GeV) [10]. 

Particle identification is performed using likelihood functions to parametrize the sharpness of the 
Cherenkov rings, which are more diffused for electrons than for muons. The zenith angle distributions 
for e-like and /i-like sub- GeV and multi- Ge V events are shown in Fig. 3 left. The MC problem described 
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Figure 3: Left: SK zenith distributions (black points) for e-like and /i-like sub-GeV, multi-GeV, throughgoing 
and stopping muons. Solid lines are no oscillation MC predictions. Right: L/E u distribution for /i-like 
events. 



before exists also in SK [10]: the e-like events were in agreement with the HKKM95 MC predictions for no- 
oscillations; they are higher than the HKKM01 non oscillated MC. For /i-like events, the new MC predictions 
are low: to reduce these problems the normalization is a free parameter. 

The ratios e-like events/MC do not depend from L/E v while /i-like events/MC show a dependence on 
L/E v consistent with an oscillation hypothesis, Fig. 3 right. The overall best fit of SK data corresponds to 
maximal mixing and Am 2 = 2.5 • 10 -3 eV 2 . 

Exotic oscillations. MACRO and SuperK data were used to search for sub-dominant oscillations due 
to Lorentz invariance violation (or violation of the equivalence principle). In the first case there could 
be mixing between flavor and velocity eigenstates. 90% c.l. limits were placed in the Lorentz invariance 
violation parameters |Av| < 6 • 10" 24 at sin 2 20 v =0 and \Av\ < 4 • 10" 26 at sin 2 26 v = ±1 [22]. Neutrino 
decay could be another exotic explanation for neutrino disappearence; no radiative decay was observed 
[23]. 

4 Long baseline neutrino beams and experiments 

Neutrino physics has opened new windows into phenomena beyond the Standard Model of particle physics. 
Long baseline neutrino experiments may allow further insight into v physics. The first long baseline v beam 
was the KEK to Kamioka (K2K) beam, the 2 nd was the Fermilab to the Soudan mine beam (NuMi). CNGS 
[24] was commissioned in 2006 and started sending neutrinos to the GS Lab. 

The CNGS beam. Fig. 4 left shows the main components of the beam at CERN [24]. A 400 
GeV beam is extracted from the SPS and is transported to the target. Secondary pions and kaons are 
focused into a parallel beam by 2 magnetic lenses, called horn and reflector. Pions and kaons decay into 
and \i in a decay pipe. The remaining hadrons are absorbed in the hadron stop. The /i's are monitored in 
2 detectors. 

Fig. 4 right shows the path of the CNGS z^'s from CERN to GS. It also shows the synchronization 
via GPS of the atomic clocks at CERN and GS. The beam is optimised for producing a maximum number 
of CC v T interactions in OPERA. Fig. 5 left shows the underground layout of the SPS and of CNGS at 
CERN; Fig. 5 right shows the scheme of the SPS operation during a test run. The mean energy is 17 
GeV, the contamination ~2%, the v e (i7 e ) < 1% and the number of v T is negligible. The muon beam size 
at the 2 nd muon detector at CERN is a ~1 m; this gives a beam size at GS of a ~1 km. The first low 
intensity test beam was sent to GS in August 2006 and 3 detectors (OPERA, LVD and Borexino) obtained 
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Figure 4: Left: The main components of the CNGS neutrino beam at CERN. Right: Sketch of the 730 km 
neutrino path from CERN to Gran Sasso and GPS time selection of events. 

their first events. The low intensity CNGS was stable and of high quality. The shared SPS beam sent a 
pulse of 2 neutrino bursts, each of 10.5 /isec duration, separated by 50 ms, every 12 s. A higher intensity 
beam expected for October 2006 did not happen because of a water leak at CERN. In September 2007 a 
2 nd test was successful, but a high intensity run was cancelled because of cooling problems. 




Figure 5: Left: Underground layout of the SPS and of the CNGS beam at CERN. Right: Scheme of the 
SPS operation at CERN during the test runs. 

The K2K experiment confirmed the results of atmospheric oscillations [13]. 

The MINOS experiment on the NuMi low energy neutrino beam is a large magnetised steel scintillator 
tracking calorimeter, complemented by a similar near detector and a calibration detector. The experiment 
obtained important results which confirm the atmospheric v oscillation picture with maximal mixing and 
Am 2 = 2.38 eV 2 [14]. 

We shall now concentrate on the experiments at GS which may use the CNGS beam. 

Borexino is an electronic detector designed to study solar z/ e 's, in particular those coming from Be 7 
decays in the center of the sun [25]. The important part of the detector is a large sphere, whose central 
part was filled with liquid scintillator at the beginning of 2007. The outer sphere was started to be filled 
with water in 2006; during the first run only few meters of water were in the sphere and only few events 
were observed. 

LVD is an array of liquid scintillators with a total mass of 1000 t; it is designed to search and study 
9 e 's from gravitational stellar collapses [16]. LVD plans to be a neutrino flux monitor of the CNGS beam. 
LVD has 3 identical "towers", each containing 8 active modules; a module has 8 counters of 1 x 1 x 1.5 m 3 , 
filled with 1.2 t of liquid scintillator. CNGS v^s are observed through the detection of muons produced 
in neutrino CC interactions in the surrounding rock or in the detector and through the detection of the 
hadrons produced in neutrino NC/CC interactions inside the detector. In the 2006 test run LVD counted 
50-100 muons per day and recorded ~500 events [16]. 

OPERA [15] is a hybrid-emulsion-electronic detector, designed to search for the v ^ < — > v T oscillations 
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in the parameter region indicated by the atmospheric neutrinos, confirmed by the K2K and MINOS 
experiments. The v T appearance will be made by direct detection of the r lepton, from v T CC interactions 
and the r lepton decay products. To observe the decays, a spatial resolution of ~ 1 fim is necessary; this 
is obtained in emulsion sheets interspersed with thin lead target plates (Emulsion Cloud Chamber (ECC)). 
OPERA may also search for the subleading v ^ < — > v e oscillations and make a variety of observations 
using its electronic detectors. The OPERA detector, Fig. 6, is made of two identical supermodules, each 




Figure 6: Layout of the OPERA detector 



consisting of a target section with 31 target planes of lead/emulsion- film modules ("bricks"), of a scintillator 
tracker detector and of a muon spectrometer. The final target mass is 1.55 kt. 

The first electronic subdetector is an anticoincidence wall to separate \i events coming from interactions 
in OPERA from those in the material and rock before OPERA. 

The target tracker is made of scintillator strips, each 7 m long and of 25x15 mm 2 cross section. A 
wavelength shifting fibre of 1 mm diameter transmits the light signals to both ends. The readout is done 
by 1000 64 channel Hamamatsu PMTs. 

The muon spectrometer consists of 2 iron magnets instrumented with Resistive Plate Chambers (RPC) 
and drift tubes. Each magnet is an 8 x 8 m 2 dipole with a field of 1.52 T in the upward direction on one 
side and in the downward direction on the other side. A magnet consists of twelve 5 cm thick iron slabs, 
alternated with RPC planes. 

The precision tracker measures the muon track coordinates in the horizontal plane. It is made of drift 
tube planes, placed in front and behind each magnet and between the 2 magnets. The muon spectrometer 
has a Ap/p < 0.25 for muon momenta < 25 GeV/c. Two 45° crossed planes of glass RPC's (XPC's) are 
installed in front of the magnets. 

The DAQ system uses a Gigabit network of 1150 nodes. To match the data of the different subdetectors 
a "time stamp" is delivered by a clock using the GPS. Also the synchronization with the beam spill is done 
via GPS. The commissioning of each electronic detector was made with CR muons and with the CNGS at 
reduced intensity. 

Nuclear emulsions and their scanning. The basic target module is a "brick", consisting of 56 lead plates 
(1 mm thick) and 57 emulsion layers. A brick has a size of 10.2 x 12.7 cm 2 , a depth of 7.5 cm (10 radiation 
lengths) and a weight of 8.3 kg. Two additional emulsion sheets, the changeable sheets (CS), are glued on 
its downstream face. The bricks are arranged in walls. Within a brick, the achieved spatial resolution is 
< 1 /im and the angular resolution is ~2 mrad. Walls of target trackers provide the v interaction trigger 
and the identification of the brick in which the interaction took place. 

The bricks are made by the Brick Assembling Machine (BAM), which consists of robots for the 
mechanical packing of the bricks. The BAM, installed in the GS lab, produces 1 brick every ~2 minutes. 
The bricks are handled by the Brick Manipulator System (BMS), made of two robots, each operating at 
one side of the detector. 
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A fast automated scanning system with a scanning speed of ~ 20 cm 2 /h per emulsion (44 /im thick) 
is needed to cope with daily analyses of many emulsions. This is a factor of 10 increase with respect to 
past systems. For this purpose were developed the European Scanning System (ESS) in Europe [26] and 
the S-UTS in Japan [27]. An emulsion is placed on a holder and the flatness is guaranteed by a vacuum 
system. By adjusting the focal plane of the objective, 16 tomographic images of each field of view are 
taken at equally spaced depths. The images are digitized, converted into a grey scale, sent to a vision 
processor and analyzed to recognize sequences of aligned grains. The 3-dimensional structure of a track 
in an emulsion layer (microtrack) is reconstructed by combining clusters belonging to images at different 
levels. Each microtrack pair is connected across a plastic base to form the base track. A set of connected 
base tracks forms a volume track. 




Figure 7: Left: display of one event with the v interacting in a "brick" in the first module, as seen by the 
electronic detectors; right: the interaction vertex seen in nuclear emulsions. 

In the 2006 test run OPERA recorded 319 events, consistent with the 300 events expected from the 
delivered integrated intensity. A fit to the 9 angle distribution of these reconstructed events on-time with 
the beam yielded a mean muon angle of 3.4° in agreement with the value of 3.3° expected for v ^ originating 
from CERN and travelling under the earth surface to the GS underground halls. During the same test run, 
a test of the CS procedure was performed, using an emulsion detector plane inserted in the SM2 target. 9 
muons produced by neutrino interactions in the rock surrounding the detector crossed the CS plane surface. 
The test proved the capability of going from the centimetre scale of the electronic tracker to the micrometric 
resolution of nuclear emulsions [15]. During the short 2007 test run several interactions in the lead bricks 
positioned in the target section of SM 1 (>40% of total mass) were detected by the electronic detectors 
and confirmed by tracks observed in the CS: the event vertexes were observed in the emulsions analysed by 
microscopes: Fig. 7 left shows an interaction in the electronic detectors, Fig. 7 right the vertex region in 
the emulsions. This test confirmed the validity of the methods to associate electronic detectors to nuclear 
emulsions. 

5 Conclusions 

The atmospheric neutrino anomaly became, in 1998, the atmospheric neutrino oscillation hypothesis with 
maximal mixing and Am^ ~ 2.4 • 10 -3 eV 2 . 

It was later confirmed with more data and by the first two long baseline experiments. All experiments 
agree on maximal mixing, while the Am^ are: 

Soudan-2 5.2 , MACRO 2.3 , SK 2.5 , K2K 2.7 , MINOS 2.38 -lO" 3 eV 2 . 

These results come from disappearance experiments (y^ — > v^). An appearance experiment (y^ < — ► v T ) 
with v T detection would solve conclusively the situation. 
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Other more exotic scenarios have been investigate like: Lorentz invariance violation, v radiative decays, 
etc. Stringent limits were established for these hypotheses. 
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